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Abstract. We review the theory of how the shapes of the X-ray isophotes 
probe both the shape and radial distribution of gravitating matter in el- 
liptical galaxies in a way that is more robust than the traditional spherical 
approach. We summarize and update the previous analyses of X-ray ob- 
servations of NGC 720 and NGC 1332 and describe preliminary results 



1. Introduction 

The soft X-ray band is particularly useful for studying the mass distributions 
of elliptical galaxies, because the emission over energies ~ 0.5 — 2 keV arises 
primarily from hot gas with temperature T ~ 10 7 K; i.e. for the more massive 
ellipticals having ratios of X-ray to blue-band luminosity log 10 L x /Lb >~ —2.7 
(Canizares, Fabbiano, & Trinchieri 1987). The small mean free path (~ ^ kpc) 
indicates that the gas is a collisional fluid with an isotropic pressure tensor. 
The short sound crossing time (<~ 10 7 yr) further indicates that, to a good 
approximation, the gas is in hydrostatic equilibrium. Any streaming motions due 
to cooling flows or supernova-driven winds should be very subsonic in massive 
ellipticals; for a recent review see Sarazin (1997). 

The traditional method for obtaining the masses of ellipticals from X-rays 
assumes spherical symmetry and that the X-ray emission arises from a single- 
phase, non-rotating, ideal gas in hydrostatic equilibrium (Fabricant, Lecar, & 
Gorenstein 1980), 



where p g is the mass density of the hot gas. The temperature profile, T(r), 
figures prominently in this equation, and has historically hindered reliable mass 
determinations for most ellipticals. Even at present accurate temperature pro- 
files exist for only a few galaxies and interpretations of these is hampered by the 
possibilities of cooling flows and multi-temperature gas (Buote & Fabian 1997). 



for NGC 3923. 




(1) 



1 



Moreover, if B fields are important then an additional term needs to added to 
equation (|l|). Clearly it is important to find a way to determine M grav in a 
manner that is insensitive to T(r) and other poorly known properties of the gas. 

2. A Geometric Test for Dark Matter 

The shapes of the X-ray isophotes probe M grav without requiring the equation 
of hydrostatic equilibrium to be solved explicitly. For a galaxy of arbitrary 
shape the geometric properties of the hydrostatic equation, Vp g = — p g X7&, 
require that the gas pressure, p g , density, and gravitational potential, have the 
same shapes in three dimensions (3-D) independent of the temperature profile 
of the gas; this can be shown by taking curls of the hydrostatic equation. (One 
replaces <3? with the appropriate effective potential, 3> e //> if the gas is rotating 
appreciably.) If the gas is also adequately described by a single-phase, p g = 
Pg(p g ,T), then T has the same 3-D shapes as p g , p g , and <£. Since the X-ray 
emissivity is a function only of these quantities, j x oc p 2 g K(T) (where A(T) is the 
intrinsic plasma emissivity), we arrive at the key property (Buote & Canizares 
1994, §3.1; 1996a, §5.1), 

X-Ray Shape Theorem The X-ray emissivity and gravitational potential have 
the same 3-D shapes independent of the temperature profile of the gas. 

This relation between the shapes of j x and $ is more robust to issues like B 
fields and cooling flows than is equation (jl]). Although in principle this relation 

may be affected by arbitrary B fields, it is unaffected if p m ag oc p g as is often 
assumed (e.g., Loeb & Mao 1994). Similarly, simple models of cooling flows with 
mass dropout (White & Sarazin 1987) just add another term to j x which is only 
a function of p g and T, and thus do not disturb the relation. A possible area 
of concern for this relation (as well as for equation |l]) is if the gas is strongly 
multi-phase (see §|5]). 

The X-ray Shape Theorem allows a robust "Geometric Test" for dark matter 
in ellipticals and galaxy clusters (Buote & Canizares 1994, 1996a). That is, the 
hypothesis that gravitating mass follows the optical light can be rigorously tested 
by comparing the shape of the potential, constructed from a constant M/L 
model with the shape of j x obtained from X-ray imaging data. The shape of 
§L depends on both the shape and concentration of pi, the mass density of the 
constant M /L model. 

This dependence on mass concentration is illustrated in Figure [l] where the 
potential exterior to a thin ellipsoidal shell with ellipticity, e p = 0.5, is shown. 
The exterior isopotentials are confocal to the shell and thus their ellipticities (e<j>) 
fall to nearly zero at only a few mean radii (r) from the shell. (The dotted line 
is e = 0.5 at 5r for comparison.) This is an example of the general property that 
the higher order potential multipoles rapidly decay with increasing distance from 
a centrally concentrated mass and give way to the spherical monopole. Hence, 
the ellipticity gradient of the X-ray isophotes, which indicates a related gradient 
in the potential, probes the radial mass distribution; see Buote & Canizares 
(1996b) for an application of this idea to clusters. 

If the constant M/L model cannot produce the X-ray isophote shapes for 
a galaxy, then one can add a dark component until agreement is achieved. Rea- 
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Figure 1. Potential shapes and mass concentration. 

sonable choices for this dark component are a halo having either a NFW density 
profile (see Navarro, these proceedings) or a Hernquist (1990) profile with the 
same ellipticity as the light. This approach allows for a lower limit on the ratio 
of dark to luminous matter, M^m/Ml, to be determined independent of T(r). 

Finally, this "Geometric Test" can test the viability of alternative grav- 
ity theories like MOND (Milgrom 1983), perhaps the most successful of its 
kind. Typically these alternative theories are devised to explain radial mani- 
festations of dark matter in galaxies (e.g., spiral galaxy rotation curves). Buote 
& Canizares (1994) showed that potential shapes in MOND are essentially iden- 
tical to those in Newtonian theory. Hence, if a constant M/L model in Newto- 
nian theory is unable to produce the observed X-ray isophote shapes of a galaxy, 
neither will the corresponding model in MOND. 

3. Observational Results for the Geometric Test 

The isolated E4 galaxy NGC 720 is the prototype for X-ray shape analysis (see 
Buote & Canizares 1994, 1996c, 1997; Romanowsky & Kochanek 1997). In 
Figure § we show the X-ray contours of the 20ks ROSAT PSPC and 57ks HRI 
observations of NGC 720 overlaid on the digitized POSS image ; 10" ~ lh^ 
kpc. (The X-ray images have been smoothed for display with Gaussians having 
half the widths of the respective PSFs; FWHM PSF is ~ 30" for PSPC and is 
~ 4" for HRI.) The X-ray isophotes are clearly flattened and the orientiations 
appear to be offset from the optical major axis at large radii (~ 100"). Here we 



concentrate on the ellipticity of the X-ray surface brightness and defer discussion 
of the orientations to §^]. 

Because of the relatively low S/N data of the X-ray images, we use a moment 
method to compute the ellipticity of the X-ray surface brightness analogous to 
computing the principal moments of inertia within an elliptical aperture (see 
Buote & Canizares 1994). The X-ray ellipticity profiles are shown in Figure ||| 
with la error bars from 1000 Monte carlo realizations of the images (see Buote 
& Canizares 1997); also displayed are the B-band isophotal ellipticities. The 
ellipticity of the PSPC data is largest and most significant for semi-major axes 
a ~ 90" — 100", while the HRI measures similar ellipticity down to a ~ 40". 

By assuming M oc Lb we compute the potential <&b] the isopotential ellip- 
ticities appear in the figure. The ellipticities of &b are considerably smaller than 
for the B-band light because Lb is highly centrally concentrated: R e ~ 50", core 
~ 4". That is, the monopole dominates &b for r >~ 10". Since the ellipticities 
of ^b fall well below those computed for the X-ray data, the M oc Lb model 
would appear to fail. However, to rigorously implement the Geometric Test we 
must deproject the X-ray data to compare to the shape of &b- 

The procedure we adopt, appropriate for the relatively low S/N X-ray data 
(~ 1500 counts), is to first jointly fit a simple model to the radial surface bright- 
ness of the PSPC and HRI data (e.g., a (5 model). The best-fit model is depro- 
jected to 3-D and assigned the ellipticities of &b- By projecting this ellipsoidal 
model back onto the sky plane (and adjusting the free parameters to maintain 
a best fit of the radial surface brightness profile), we obtain a M oc Lb model of 
the X-ray surface brightness. As above, we perform 1000 Monte carlo realiza- 
tions of this model and compute moment ellipticities, e x ($s), analogous to the 
data to arrive at the la error bars in the bottom panels of Figure [2|. 

The observed PSPC and HRI ellipticities exceed the 6 x (&b) generated by 
the M oc Lb model at greater than the 3a level, thus indicating the need 
for flattened dark matter. Moreover, even if Lb is allowed to have ellipticity 
~ 0.6, consistent with the flattest ellipticals observed and expected from stability 
arguments, this discrepancy cannot be entirely resolved. Hence, in NGC 720 the 
dark matter must be both flattened and more extended than Lb, a conclusion 
independent of the gas temperature profile (Buote & Canizares 1994, 1997). 

We have also applied the Geometric Test to PSPC data of the fairly isolated 
E7/S0 galaxy NGC 1332 (left panels of Figure |) and to PSPC data of the classic 
shell galaxy - the isolated E4, NGC 3923 (right panels of Figure ||) . Both of these 
galaxies have lower S/N images than NGC 720 and have smaller measured X-ray 
ellipticities. Similar to NGC 720, the quite flattened optical light cannot produce 
flat enough X-ray isophotes. For NGC 1332, the observed PSPC ellipticities 
exceed the e x (&i) generated by the M oc Lj model at the 90% level and, similar 
to NGC 720, both flattened and extended dark matter are indicated (Buote & 
Canizares 1996a). The M oc Lb model for NGC 3923 disagrees only at the 80% 
level. Although this latter case is marginal, the character of the discrepancy is 
the same (i.e. flattened and extended dark matter needed - Buote &: Canizares 
1997, to be submitted to ApJ). 
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Figure 2. NGC 720 PSPC (L) and HRI (R) 
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Figure 3. NGC 1332 (L) and NGC 3923 (R) 




X-Ray and Optical Ellipticities 





50 100 50 100 

a (arcsec) a (arcsec) 

Geometric Test for DM 



1 1 1 1 1 

; Boxes - 


1 

-> 


1 1 1 

PSPC 


i 


1 — 


Error 


Bars 


— > : 




D □ 






□ 




: 






' [ 


1 " 




i 






1 






3 


50 






100 







— i 1 1 1 1 1 1 1 1 1 1 1 r— 

Boxes --> PSPC 
Error Bars — > e x (^ R ) 



a (arcsec) 



50 100 

a (arcsec) 



4. Detailed Hydrostatic Models 



The Geometric Test robustly examines the need for dark matter and provides 
some constraints on its shape and extent. To obtain detailed information on the 
range of allowed shapes and radial profiles for the dark matter the hydrostatic 
equation must be solved explicitly. Following the pioneering approach of Binney 
& Strimple (1978), one may solve the hydrostatic equation for the gas density 
by assuming a single-phase, non-rotating ideal gas, 

p g = f- l exp^-rjf- 1 V$-dzj, (2) 

where V = p,m p <&(Q)/kBT(0) and where, e.g., p g = p g (x)/ p g (0). Starting with a 
model for the mass distribution one computes and then, in conjunction with 
a measurement (or assumption) for the temperature profile, one computes p g . 
The model X-ray surface brightness is then obtained by integrating p^ along the 
line-of-sight. This is appropriate since A(T) convolved with the ROSAT spectral 
response is effectively constant for the relevant range of temperatures. 

Although not immediately apparent from equation (Q), constraints on the 
shape of the mass distribution using this method are quite insensitive to tem- 
perature gradients (though not independent of them like the Geometric Test). 
Strimple & Binney (1979) were the first to show that isothermal and adiabatic 

(p g oc p g ^ 3 ) temperature profiles lead to approximately similar X-ray isophote 
shapes. Hence, an isothermal gas may be reasonably assumed if one is only 
interested in measuring the shape of the gravitating mass. 

The results for the ellipticity of the gravitating matter in NGC 720, NGC 
1332, and NGC 3923 are shown in Figure H]. These correspond to oblate 
spheroidal models with mass density, p ~ r _ , and for an isothermal gas (see 
Buote & Canizares 1994, 1996a, 1997). (The results do not vary much for other 
models.) The intensity weighted optical ellipticity (e.g., < B >) and the maxi- 
mum value (e.g., B max ) are also shown. One can see that the ellipticity of the 
gravitating matter is at least as large as that of the light in all three cases, with 
some indication (especially NGC 720) that the dark matter is more elongated. 

Finally, the scale lengths of the gravitating matter are ~ 5 times that of 
the stellar luminosity for NGC 720, and ~ 3 times for NGC 1332 and ~ 2 
times for NGC 3923. These scale length ratios are quoted for Hernquist models 
to insure a consistent comparison between the light and mass, and the precise 
values depend on the temperature profile. 

5. Caveats 

The principal caveats associated with the analysis of X-ray isophote shapes 
described above are, (1) multi-phase gas, (2) discrete sources, (3) rotation, and 
(4) environmental effects. There is new evidence from ASCA that the hot gas in 
the brightest ellipticals consists of at least two temperatures (Buote & Fabian 
1997). The single-phase X-ray shape analysis is insensitive to multiple gas phases 
if either (a) a single phase dominates the emission, (b) (see j p^| ) the emission is 
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Figure 4. Ellipticity of Gravitating Mass 



adequately described as a simple superposition of an ambient phase and a mass 
dropout term (White & Sarazin 1987), or (c) the phases have similar spatial 
distributions so that the emission-weighted (or mass-weighted) gas density and 
temperature are good descriptions of the emissivity; i.e. j x oc (p g ) 2 A((T)). 
These mean quantities appear to provide a good description of general multi- 
phase cooling flows in clusters (Thomas, Fabian, & Nulsen 1987). 

Discrete sources should contribute to the soft X-ray emission of early-type 
galaxies. In fact, the ASCA spectra of NGC 720, NGC 1332, and NGC 3923 
do require a second thermal component with T >~ 5 keV suggestive of discrete 
sources (Buote & Canizares 1997; Buote & Fabian 1997). However, the large 
inferred temperatures of these hard components in these relatively low S/N 
galaxies may be a fitting artifact since the highest S/N ellipticals have "hard" 
components with T <~ 2 keV which instead indicate another phase of hot 
gas (Buote &; Fabian 1997). Moreover, the spatial distribution of the discrete 
component has not been accurately measured. However, assuming the discrete 
sources contribute to 20% of the ROSAT emission of NGC 720 and that they 
are distributed like the optical light, Buote & Canizares (1997) have shown that 
the inferred ellipticity of the gravitating mass decreases by only ~ 0.05. 

It is unlikely that rotation plays an important role in the dynamics of NGC 
720 and NGC 3923 since the stellar rotation is negligible. Theoretically, even 
without strong stellar rotation one may (Kley & Mathews 1994) or may not 
(Nulsen, Stewart, & Fabian 1984) expect a rotating cooling flow to develop de- 
pending on whether angular momentum of the gas is conserved. Highly flattened 
X-ray isophotes indicative of a cooling disk have not been observed. Finally, it 
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Figure 5. X-Ray Position Angle Twist in NGC 720 



is unlikely that environmental effects substantially distort the X-ray isophotes 
in the isolated ellipticals studied since for r <~ 100" no large asymmetries 
indicative of ram pressure or tidal distortions are observed. 



6. X-ray position angle twist in NGC 720. 

The major axes of the X-ray isophotes of NGC 720 (Figure ||) are offset from the 
optical major axes. In the left panel of Figure |H| we show the position angle (PA) 
profile of the HRI data with the PAs of the PSPC and optical data indicated 
for comparison (Buote & Canizares 1996c). For r <~ 60" the X-ray isophotes 
are aligned with the optical and then twist away to a maximum misalignment 
of ~ 30 degrees at r ~ 100". This misalignment provides further geometric 
evidence for dark matter. Buote & Canizares (1996c) attempted to explain this 
PA twist with a simple triaxial halo model where the axial ratio varied as a 
power law in radius but were unable to simultaneously produce the large twist 
while maintaining the relatively large X-ray ellipticities. 

Romanowsky & Kochanek (1997) have recently examined more detailed 
physical triaxial models to try to explain the NGC 720 X-ray data. Although 
the best model of Romanowsky & Kochanek (solid line in right panel of Figure 
||) is able to produce a large twist, the twist occurs at very small radii and is 
thus unable to explain the abrupt twist at ~ 60" seen in the HRI data. An 
intrinsic misalignment of the stars and dark matter halo may be implicated. 



7. Conclusions 



X-ray isophote shapes probe both the shape and concentration of gravitating 
mass essentially independent of the temperature profile of the gas. (The Ge- 
ometric Test is completely independent of T(r).) For the early-type galaxies 
studied (NGC 720, NGC 1332, NGC 3923), dark matter is required to explain 
the elongated isophotes. (MOND is unable to escape this manifestation of dark 
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matter.) The dark matter is at least as flattened as the optical light and is also 
more extended for each of the galaxies. The X-ray position-angle twist in NGC 
720 appears to suggest an intrinsic misalignment of stars and dark matter halo 
rather than a strongly triaxial system. 

This X-ray shape analysis is applicable to isolated early-type galaxies with 
typical log 10 L x /Lb >~ —2.7 (i.e. X-rays from hot gas). The next generation X- 
ray satellites soon to be flown, AXAF and XMM, will have the spatial resolution 
and effective area to precisely map X-ray isophote shapes for a statistically large 
sample of galaxies. X-rays are the most promising means on the horizon for 
obtaining a large sample of gravitating mass shapes and radial mass profiles of 
early-type galaxies which should significantly enhance our understanding of the 
structure and formation of these systems (e.g., Sackett 1996; de Zeeuw 1997). 
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